Frontiers in Neuroimaging

Multisite functional Magnetic Resonance Imaging (fMRI) studies are rapidly becoming the norm to allow the acquisition of large samples required to perform advanced statistical techniques, e.g. machine learning. However, acquiring data at multiple sites includes methodological and technical challenges due to the difference in the facilities of each site, such as scanner manufacturer and model. These new challenges add to the already well-known challenges of fMRI data, including scanner-related noise, participant movement, and physiological confounds. To ensure the optimal quality of data, a balance needs to be achieved between standardizing data preprocessing across sites and optimizing within-site data preprocessing. To define the optimal preprocessing pipeline including standard steps and additional denoising technique for our multisite dataset, we first tested 3 commonly used preprocessing pipelines, i.e. fMRIPrep, FSL, and CONN. These pipelines all include standard steps, e.g. motion correction, temporal filter, etc. In a second step, to further improve signal quality, we tested the efficiency of 3 additional denoising techniques on the output of the data preprocessed with the previously defined pipeline. These techniques included aCompCor, FSL FIX and ICA-AROMA. Signal quality was quantified as temporal signal-to-noise ratio (tSNR) across and within sites. Our results show that the performance of the preprocessing pipeline and denoising technique varies between sites. FSL yielded the highest tSNR at two sites and fMRIPrep at the third, revealing a significant site x pipeline interaction. An FSL-based pipeline with minimal adjustment to accommodate site-specific parameters and followed by denoising using a single FSL FIX classifier, which was custom-trained across sites, achieved the highest quality of signal in our data and yielded the greatest consistency in improved data quality across sites. Because of the potential of this pipeline to adequately identify and address site-specific noise, while remaining constant across sites, we selected it as optimal preprocessing pipeline for our dataset.

Background: MRI has increasingly been explored as a biomarker for detecting structural and functional brain changes. For clinical decision-making, it is crucial to validate observed changes in MRI indices at the individual level. The uncertainty in longitudinal MRI indices can be quantified using the repeatability coefficient (RC). Methods: Twenty healthy controls (10 males, 10 females) underwent two test-retest sessions of structural magnetic resonance imaging (MRI) and resting-state functional MRI (rs-fMRI) on the same day, separated by a 30-minute interval. RC values and their 95% confidence intervals (CI) were estimated for subcortical volumes, cortical thickness, and within-network functional connectivity. Additionally, 33 patients with mental health disorders underwent MRI before and after 20 sessions of transcranial magnetic stimulation (TMS). Percentage changes in MRI-derived indices were assessed at the individual level, with changes exceeding the RC threshold considered indicative of true change beyond measurement uncertainty. Results: The RC showed measurement variability in subcortical volumetric in the range of 8% to 17.5% for caudate and left amygdala, respectively. For cortical thickness, the RC was measured between 3.5% and 16.5% for the left occipital pole and the left temporal pole, respectively. The RC% for fractional anisotropy (FA) measures were variable between 11.3% (the left middle cingulum) and 62.9% (the right anterior cingulum). For within-network connectivity, the RC was measured in a range of 9.7% and 29.4% for sensorimotor and visual networks, respectively. TMS-treated patients exhibited no changes beyond the RC in almost all subcortical volumes and within-network connectivity. The most frequent changes beyond the uncertainty were observed in FA measures, particularly in the posterior cingulum, where 17 out of 23 patients exhibited clinically meaningful alterations. Conclusion: Structural brain features extracted from MRI demonstrated high reliability. Among all measures, FA, reflecting white matter integrity, was most sensitive in detecting neural changes following TMS, highlighting its potential utility as a treatment-responsive biomarker.

1.IntroductionUnderstanding mid-term test-retest reliability and within-subject variability is important for interpreting changes observed in longitudinal and intervention studies. The reliability of resting-state functional magnetic resonance imaging (rs-fMRI) is known to vary across measures and brain regions. However, how reliability differs across functional networks and connectivity-and amplitude-based measures, and whether multi-echo acquisition and processing modify these patterns, remain incompletely characterized. MethodsTwenty-two healthy volunteers underwent two rs-fMRI sessions 15.7 {+/-} 4.0 days apart on a Compact 3T scanner. Multi-echo, middle-echo, and independently acquired single-echo datasets were compared, with multi-echo independent component analysis additionally evaluated as a denoising approach. Functional connectivity (FC) and three amplitude-based measures were evaluated using the Schaefer 400 parcellation. Reliability was systematically assessed using intraclass correlation coefficient (ICC), within-subject standard deviation (wSD), and systematic bias at edge or regional, and network levels. ResultsAcquisition-dependent differences in reliability were generally modest. Multi-echo acquisition and processing increased functional connectivity strength and the magnitude of amplitude-based measures and improved inferior cortical coverage, but these enhancements did not consistently translate into substantially higher ICC or lower wSD. In contrast, reliability showed clear network-dependent differences. FC reliability varied markedly across network pairs and was not explained by connectivity strength alone; pairs involving the default mode and control networks generally showed more favorable profiles than several somatomotor and visual network pairs. Fractional amplitude of low-frequency fluctuations (fALFF) also showed network-dependent reliability, with the most favorable regional reproducibility observed in the default mode and control networks and lower reproducibility in the somatomotor and visual networks. ConclusionThese findings provide practical mid-term reliability benchmarks for rs-fMRI on a Compact 3T scanner and show that measurement stability varies more clearly across measures and functional networks than across acquisition approaches. Key pointsO_LIMid-term test-retest reliability varied more clearly across resting-state measures and functional networks than across acquisition and processing approaches. C_LIO_LIMulti-echo acquisition and processing enhanced functional connectivity strength, amplitude-based signal magnitude, and inferior cortical coverage but did not consistently improve reliability. C_LIO_LIFunctional connectivity strength and fractional amplitude of low-frequency fluctuations showed distinct network-specific reliability profiles, with more favorable reproducibility in default mode and control networks than in several somatomotor and visual networks. C_LI

Despite the growing use of functional near-infrared spectroscopy (fNIRS) hyperscanning to record brain activity simultaneously from interacting individuals in naturalistic settings, most analyses quantify functional connectivity separately for each channel pair. The resulting collection of pairwise estimates is difficult to integrate into a network-level characterization of intra- and inter-brain organization. Here, we present an open, configuration-driven Python toolkit that transforms preprocessed fNIRS hyperscanning time series into functional connectivity graphs. The toolkit constructs a bipartite inter-brain network for each dyad and separate intra-brain networks for each participant, computes node- and graph-level measures, and exports adjacency matrices, edge lists, analysis-ready summary tables, reproducibility metadata, and standardized visualizations. Dataset-specific parameters, including directory structure, participant naming, channel selection, epoch extraction, and edge-retention criteria, are defined in a human-readable YAML configuration file, enabling the same workflow to accommodate differently organized datasets without changes to the source code. We illustrate the pipeline using a representative recording from a mother-infant fNIRS hyperscanning dataset and present the resulting network outputs. The toolkit provides a reproducible framework for moving from pairwise functional connectivity estimates to network-level analyses of dyadic and individual brain organization.

Background. Cerebrovascular reactivity (CVR) is a sensitive marker of vascular health, most commonly measured with BOLD fMRI during hypercapnic gas challenges. Hypoxia is a complementary but under-characterized vasoactive stimulus. We characterize a combined hypercapnic-hypoxic gas-challenge paradigm at 7 Tesla BOLD fMRI and its cross-session reproducibility. Methods. Eleven healthy adults were scanned twice at 7 Tesla during an 18-minute multi-echo, multiband BOLD acquisition, receiving fixed-inspired hypercapnic (5% CO2) and hypoxic (10% O2) gas boluses separated by medical-air recovery periods, with end-tidal gas and peripheral saturation monitoring. Hypercapnic and hypoxic CVR were estimated using three trace-based general linear models: end-tidal gas, peripheral saturation, and Severinghaus-derived arterial saturation; and a finite impulse response model. Cross-session reproducibility was assessed by region-of-interest intraclass correlation. Results. The paradigm produced robust, condition-specific gray-matter BOLD responses that reproduced across sessions. The Severinghaus-derived saturation regressor generalized significantly better than end-tidal and peripheral regressors under cross-session validation and outperformed the finite impulse response benchmark in 58% of gray-matter voxels. Reproducibility across cortical, subcortical, and brainstem atlases was moderate-to-high (median cortical intraclass correlation 0.71). Conclusion. A combined hypercapnic-hypoxic paradigm yields reproducible voxelwise cerebrovascular reactivity maps at 7 Tesla; Severinghaus-derived arterial saturation is the most generalizable regressor for hypoxic reactivity.

Multi-echo functional magnetic resonance imaging (fMRI) acquires signals with distinct contrast profiles across echo times (TEs). Empirical evidence suggests that different TEs capture distinct signal contributions with varying sensitivity to BOLD and non-BOLD processes from long to short TEs (Stroman, 2002; Krishnamurthy, 2023; Dong, 2024). Previous literature has not yet examined how brain states extracted from temporal dynamic functional connectivity (dFNC) evolve across TEs. Here we present a methodological and empirical study of TE-dependent evolution of dynamic connectivity states, with a phenomenological two-component model for interpretation. Multiple echoes were acquired with an echo planar time-resolved imaging (EPTI) sequence. TE-specific components were extracted via atlas-based group ICA (Neuromark 1.0), and brain connectivity states were generated via windowed dFNC and spatial ICA clustering of windows. Within-subject and group-level states evolved across TEs in network organization and variance. States exhibited TE-dependent variance profiles, with distinct peak echo times; and group-level points of state transition were independently mapped to TE-dependent transitions in variance. We modified a bi-exponential model to simulate and model changes in the dFNC time-course across TEs. The bi-exponential model revealed high consistency with subject states (R2 fit=0.86, res. error=0.55, 86% model convergence rate), suggesting a two-component signal system contributing to dynamic TE-dependent connectivity changes in state. A significant correlation between relative model terms was estimated on a cell-wise basis across states (r(24) = -0.43, p=0.03, 95% CI [-0.70, -0.04]). This suggests that specific bi-exponential signal parameters (i.e., short FNC1,TE=0, long FNC2,TE=0, and long FNC2,decay) were associated with the observed TE-dependent variations in connectivity states across TEs. In sum, our analyses provide downstream implications for TE-aware state analysis dependent on connectivity magnitude (FNCTE=0) and rate of connectivity decay (FNCdecay), in a two-component signal model, with future implications for TE-dependent state biomarkers.

Quantitative magnetic resonance imaging (qMRI) supports both basic neuroscience and clinical research, providing the foundation for longitudinal studies of brain development, aging, and neurological disease. However, its validity rests on a little-studied assumption: that physical relocation of imaging hardware leaves underlying neuroimaging biomarkers intact. When research facilities expand or modernize, an entire scanner undergoes radical disruption--field ramp-down, physical transport, shielding reconstruction, and re-shimming. Whether qMRI metrics remain stable following such invasive procedures without introducing systematic biases is largely unquantified. Here, we systematically evaluated the impact of relocating a 3T clinical scanner over a distance of 30 km on a broad, multimodal suite of neuroimaging markers in healthy volunteers. To disentangle variance components, we benchmarked within- and cross-site stability using percent variation (PV; measurement precision), intraclass correlation coefficients (ICC; absolute reliability), and Spearman rank correlations (R2; spatial and order agreement). We show that scanner relocation introduces no greater variance than baseline within-site assessments. High-resolution brain structural morphometry from different sequences (multi-echo MPRAGE, compressed sensing MPR2RAGE) and quantitative relaxometry (compressed sensing MP2RAGE T1 mapping) demonstrated excellent stability (PV<1%, ICC=0.89-0.99, R2>0.93). Diffusion tensor (mean diffusivity: MD, fractional anisotropy: FA) and mean signal kurtosis (MSK) metrics from multi-shell diffusion MRI maintained low variability across white and deep gray matter (global cross-site ICC>0.90; global PV: MD<1.7%, MSK<1.2%, FA<3.8%; global R2: MSK>0.93; MD>0.90; FA>0.77), confirming relocation robustness for diffusion-derived microstructural markers. Resting-state functional networks displayed high cross-site reproducibility (ICC>0.90), while task-based fMRI activation peaks retained acceptable spatial and magnitude agreement (ICC=0.30-0.60) despite intrinsic physiological variability. In the absence of software or hardware upgrades, physical scanner relocation leaves qMRI biomarkers remarkably intact--providing the empirical foundation required to preserve data integrity across major facility transitions and multi-site longitudinal cohorts. Key points that summarize the paperO_LIHardware transitions do not compromise qMRI data integrity: Physical relocation of a 3T scanner (including field ramp-down, 30 km transport, and re-shimming) introduces no greater variance in quantitative MRI metrics than baseline within-site test-retest assessments. C_LIO_LIMultimodal biomarker robustness: Brain morphometry, T1 relaxometry, diffusion tensor and kurtosis, resting-state functional connectivity display comparable within- and cross-site stability reflected by different metrics (percent variation, intraclass correlation coefficient and Spearman correlation) C_LIO_LIEmpirical backing for longitudinal research: In the absence of hardware or software upgrades, multimodal qMRI biomarkers remain remarkably stable, reassuring researchers that long-term cohort studies and multi-site harmonization can survive major facility relocations. C_LI Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=175 SRC="FIGDIR/small/26363342v1_ufig1.gif" ALT="Figure 1"> View larger version (57K): org.highwire.dtl.DTLVardef@773d03org.highwire.dtl.DTLVardef@46a09borg.highwire.dtl.DTLVardef@da5bc0org.highwire.dtl.DTLVardef@1bdbf9d_HPS_FORMAT_FIGEXP M_FIG This study evaluates the impact of physically relocating a 3T MRI scanner on neuroimaging biomarkers. By comparing before and after relocation sessions on healthy volunteers, we conclude that the procedure does not introduce significant variability. Parameters for morphometry, T1 relaxometry, diffusion, and resting-state functional connectivity show high stability. Task-based fMRI showed greater within-site variability, but no increase in variability following relocation, consistent with the inherently higher variability of task-evoked responses.. The results confirm the integrity of qMRI data after infrastructure relocation, ensuring reliability for longitudinal studies. C_FIG

Background: Ultra-low field (ULF) magnetic resonance imaging (MRI) offers accessible point-of-care neuroimaging, yet its capability to support multivariate network analyses remains understudied. Methods: We systematically evaluated the feasibility and reliability of Morphometric Similarity Networks (MSN) and Morphometric Inverse Divergence (MIND) networks derived from paired high-field (HF, 3T/1.5T) and ULF (64mT) scans in a primary cohort(N=23) and an independent replication cohort (N=19). We tested the impact of super-resolution processing, cortical parcellation granularity (Desikan-Killiany vs. Destrieux), and mesh resolution (fsnative vs. fsaverage5) on edge-level cross-scanner agreement (Lin's concordance correlation coefficient, CCC) and subject-level identifiability. Results: MIND networks demonstrated higher edge-level HF-ULF agreement (mean CCC = 0.95) and ULF reproducibility (mean CCC = 0.98) compared to MSNs (mean CCC = 0.86-0.87 and 0.91-0.94, respectively). In contrast, MSNs yielded higher subject-level differential identifiability across field strengths due to a broader distinction between within- and between-subject similarity. Reconstructing MIND networks on native high-resolution meshes (fsnative) maximised cross-scanner agreement, whereas coarser meshes and finer cortical parcellations enhanced individual network identifiability. Conclusions: Structural similarity networks can be reliably constructed from ULF MRI data. Optimal analytic choices depend on study objectives: MIND on native meshes maximises cross-scanner stability, whereas MSNs excel at preserving individual uniqueness, paving the way for globally scalable connectomic research.

BackgroundSingle-voxel magnetic resonance spectroscopy (MRS) is a non-invasive method for measuring clinically and cognitively relevant metabolites. Reliable measurements require precise voxel placement across sessions and participants. We developed a scanner-console-based approach to improve voxel placement precision. MethodsIn a crossover design (n=7; six sessions each), we compared test-retest reliability of three voxel placement methods in a reference benchmark (left parietal cortex) and a technically challenging region (left ventromedial prefrontal cortex). Methods included (1) conventional anatomy-based placement, (2) mask-guided real-time positioning (MGRP), and (3) semiautomated session-locked voxel repositioning (SSVR). Resting-state MRS data were acquired using PRESS and MEGA-PRESS. Within-subject reliability of voxel placement and metabolite concentrations, namely, total N-acetylaspartate (tNAA), total Creatine (tCr), GABA (gamma-aminobutyric acid), and Glx (glutamate + glutamine) are reported using the coefficient of variation (CV), the intraclass correlation coefficient (ICC), minimal detectable change (MDC), and the spatial overlap. ResultsSSVR markedly improved voxel placement reliability, increasing spatial overlap (up to 88%) and achieving near-perfect geometric reproducibility (ICC = 0.99) compared to conventional anatomy-based placement and MGRP. SSVR improved tissue composition consistency and reduced metabolite variability in the technically challenging region (variability reduction of [~]70% tCr, [~]59% tNAA, and [~]51% Glx) while further refining already stable measurements in the benchmark region (tNAA from [~]15% to [~]10%). ConclusionBoth MGRP and SSVR improved voxel placement and metabolite measurement reproducibility compared with conventional anatomy-based placement. SSVR further enhanced within-subject reproducibility across repeated sessions, particularly in the technically challenging region, providing a robust approach for longitudinal single-voxel MRS studies.

The eigenvalue structure of resting-state fMRI (RS-fMRI) signals provides a compact representation of its variance-covariance structure, yet the information it encodes remains unclear. In this study, we introduce an eigenvalue-based framework to characterize RS-fMRI data using two features derived from the eigenspectrum: the log10-transformed first eigenvalue (log10({lambda}1)) and the slope of the log10-transformed spectrum ({beta}). We systematically evaluated these features across 14 denoising strategies using two independent datasets (China167 (83 males, 84 females; mean age +/- SD: 41 +/- 14 years old) and HCP1200 (425 males, 501 females; mean age +/- SD: 29 +/- 4 years old)). Specifically, we used singular value decomposition to obtain the eigenspectrum of the cortical BOLD signals after preprocessing and nuisance regression, and a linear model to parameterize it. We then examined the relationships between the eigenvalue parameters, head motion, and functional connectivity metrics across subjects and denoising strategies. Our key findings include: (1) log10({lambda}1) and {beta} strongly covaried with one another across subjects, denoising strategies, and datasets, indicating that they capture highly coherent aspects of the eigenspectrum; (2) both features were systematically influenced by denoising strategies; (3) within each denoising strategy, participants with greater log10({lambda}1) had greater mean framewise displacements (mFD), demonstrating sensitivity to residual motion effects, for all strategies in HCP1200 and 12/14 in China167; (4) across denoising strategies, mean log10({lambda}1) reflected the proportion of functional connectivity edges significantly associated with motion, indicating that higher log10({lambda}1) reflects more widespread motion-related contamination across large-scale functional networks; and (5) global signal regression consistently reduced log10({lambda}1), whereas spike regression had dataset-dependent effects. Together, these results establish eigenvalue signatures as robust and sensitive metrics for quantifying residual motion effects in RS-fMRI and provide a framework for evaluating denoising performance based on eigenvalue structure. The full procedure is implemented in R, and the corresponding script is available at: https://github.com/lejianhuang/EigenvalueSignature.

Magnetic resonance spectroscopic imaging (MRSI) enables non-invasive mapping of neu-rometabolites levels across the human brain. Although spectral fitting and metabolite quan-tification are increasingly supported by mature software tools, the downstream processing of quantified metabolite maps remains heterogeneous across laboratories. Here, we introduce MRSIPrep, an open-source, modular, and reproducible post-quantification framework for whole-brain MRSI. MRSIPrep standardizes quality control, tissue correction, spatial nor- malization, atlas projection, and derivative generation from quantified metabolite maps and associated quality metrics. The framework produces voxelwise, regional, and connectomics-ready outputs together with automated quality-control reports. We describe the architecture of MRSIPrep and demonstrate its utility for reproducible MRSI analysis across datasets,acquisition protocols, and downstream applications.

Purpose: Magnetic resonance spectroscopic imaging (MRSI) is a versatile technique to investigate the spatial distribution of in vivo metabolism. However, processing MRSI data is demanding, and only a few software packages support end-to-end analysis. The goal of this study was to implement fully automated, end-to-end MRSI analysis into the open-source 'Osprey-MRSI' software package. Methods: MRSI-specific analysis and visualization capabilities were implemented, building on the existing Osprey workflow. Modifications included spatial transformation and filtering operations, automated brain masking and tissue segmentation of the MRSI data, improved lipid filtering, rapid integral maps, linear-combination modeling with explicit B0 frequency-shift correction, and generation of quality-control maps and metabolic images. A fully interactive GUI and semi-interactive HTML reports provide a user-friendly way to inspect each step of the analysis. All analysis derivatives are also exported in NIfTI and NIfTI-MRS format for easy visualization and synergies with other toolboxes and modalities. Results: The automated MRSI workflow was successfully used to analyze short- and medium-TE 3T in vivo MRSI datasets from all major vendors (Philips, GE, Siemens) across multiple sites. Correct coregistration of MRSI data and MR images was validated using phantom data from each vendor and existing MRSI processing tools. Conclusion: Osprey-MRSI offers state-of-the-art methods with minimal user interaction available for non-expert users. The modularity of the workflow and the modeling algorithm will foster innovation and development of novel MRSI-specific analysis methods.

The choroid plexus (ChP) is a highly vascularized structure located in the brain ventricles that plays important roles in multiple central nervous system functions. ChP enlargement has been implicated in many neuropsychiatric conditions such as multiple sclerosis, depression, schizophrenia, and Alzheimer's disease. However, limitation in the accuracy of the current segmentation tools highlight the need for more reliable methods to quantify ChP volume in healthy and clinical human populations. The shape and size of the ChP vary drastically from person to person, and currently there are no large ChP datasets which allow comprehensive training of segmentation models. Here we introduce DataPlex, a collection of 700 structural MRI brain images where the lateral ventricle ChP has been manually labeled and checked by expert human raters. DataPlex is balanced and representative across a wide range of scanning parameters, demographics, and clinical profiles, and subject-level information will be available for each labeled image. DataPlex will be open source to support the creation of reliable segmentation tools for the human lateral ventricle ChP.

Resting-state functional MRI (rs-fMRI) studies in schizophrenia commonly rely on normalization to volumetric templates and fixed atlas parcellations derived from neurotypical populations. While these approaches enable group-level comparisons, they may obscure individual variation in cortical organization and intrinsic brain dynamics. In this study, we compared four preprocessing and parcellation strategies across two independent schizophrenia cohorts (MRI site 1, n=159; MRI site 2, n=255) to evaluate how analytic choices affect static functional connectivity and dynamic quasi-periodic pattern (QPP) measures, including default mode-dorsal attention network opposition, QPP component rank, explained variance, event rate, and associations with PANSS symptom severity. Across datasets, individualized surface-based parcellation (IndiPar) consistently detected more pronounced QPP dynamics, stronger default mode / dorsal attention network opposition, and greater explained variance of QPPs relative to atlas-based pipelines. IndiPar also produced larger and more reproducible patient-control differences in functional connectivity and QPP event-rate measures, suggesting improved sensitivity through preservation of subject-specific organization. IndiPar additionally detected a significantly increased QPP event rate and more symptom associations in patients in the larger dataset. However, associations between fMRI measures and symptom severity showed limited stability across cohorts. These findings extend previous reports of altered resting-state activity in schizophrenia, and demonstrate that preprocessing and parcellation choices substantially influence both static and dynamic rs-fMRI results. Individualized surface-based parcellation appears to better preserve subject-specific variability and improves detection of intrinsic brain dynamics. At the same time, the limited cross-dataset replication of symptom associations highlights the challenges of deriving stable brain-symptom relationships from heterogeneous psychiatric cohorts.

There is increasing demand for quantitative measurements in neuroimaging research, including absolute quantification of regional myelin density. Fast whole-brain macromolecular proton fraction (MPF) imaging enables acquisition and reconstruction of quantitative maps of MPF, a biomarker of myelin, and has been used in many research applications. To date, no standard protocols have been published for parcellating MPF maps. FreeSurfer is a widely used software package for automated brain parcellation. It is designed for use with high-resolution T1-weighted anatomical images, typically acquired with an MPRAGE sequence. MPF maps inherently provide high contrast between white and gray matter and can potentially be used as input images for parcellation, enabling simultaneous quantitative tissue characterization and volumetric assessment. However, it is unclear how the use of MPF maps as input images with FreeSurfer software affect reliability and repeatability of regional MPF and volume estimates. We evaluated three FreeSurfer parcellation workflows for this purpose. The workflows differed in the amount of neuroimaging data required and computational intensity. MPF and MPRAGE data were acquired in two separate sessions for 11 adults. Reliability was evaluated using mean relative differences in estimates, Bland-Altman analysis, and intra-class correlation coefficients (ICCs). Repeatability was assessed using within-subject coefficients of variation (CVws). All three workflows produced similarly high within-subject repeatability across scans. However, reliability in regional MPF estimates was lower in gray matter than in white matter for workflows that used MPF maps as input. Workflows that use MPF maps as FreeSurfer input may be adequate for estimating mean MPF in GM and WM of cortical parcels and in subcortical regions. However, the MPRAGE-based workflow is recommended when reliable individual-level estimates or volume estimates are required.

Transcranial magnetic stimulation (TMS) is widely used in both clinical and research settings to study and treat neurological and neuropsychiatric disorders, yet its underlying neural mechanisms remain unclear; particularly how stimulation influences both local and distant regional activity in relation to behavior. In this study, we combined single-pulse TMS with concurrent whole-head functional near-infrared spectroscopy (fNIRS) to examine hemodynamic and behavioral responses during a working memory task, with a focus on behavioral variability. Single TMS pulses were delivered to the left dorsolateral prefrontal cortex (DLPFC) while healthy participants rested; additionally, single pulses were delivered online to the left DLPFC while participants performed the working memory task. Across all participants, we observed a reliable load-dependent increase in hemodynamic activity associated with task. However, behavioral responses to TMS varied during the task. When participants were stratified into subgroups based on performance, a distinct topographic pattern emerged. During the task, TMS systematically modulated hemodynamic responses in regions including DLPFC, superior medial gyrus, precuneus, and parietal lobule, which are areas belonging to the default mode network. Moreover, the hemodynamic response during the single pulse alone sessions without any task was also found to be associated with the behavioral responses in a coherent pattern involving DLPFC, precuneus and parietal lobules. These findings suggest that variable behavioral outcomes during online TMS task are linked to distinct hemodynamic responses in a topographic pattern of local and distant regional areas.

Background and purpose: Accurate identification of cardiac MRI volume orientations is essential for reliable image interpretation and for enabling downstream automated analysis pipelines. However, orientation labels and commonly assigned manually, making the process time-consuming and prone to variability. Recent advances in vision-language models offers new opportunities for automated label generation by jointly modeling visual and semantic information. Methods: In this study, we systematically compare a conventional end-to-end transformer-based captioning model with a CLIPCap architecture for automated cardiac MRI view labelling. Both models were evaluated on a dataset of 1184 de-identified cardiac MR images spanning across six clinically relevant orientations: short axis, two-, three-, four-chamber views, left and right ventricular outflow tract views. The caption transformer was trained from scratch, whereas CLIPCap leveraged pretrained CLIP image embeddings combined with lightweight mapping network and a frozen language model. Results: Experimental results demonstrate a substantial performance advantage from CLIPCap, which achieved an overall accuracy of 98% compared to 68% for the caption transformer. CLIPCap consistently delivered high precision and recall across all orientation classes, including those of high clinical importance, while the end to end transformer exhibited unstable performance and complete failure in certain views. Conclusion: These findings highlight the benefits of leveraging large scale multimodal pretrained model for medical image labeling tasks, particularly in data limited settings. The study suggests that CLIPCap provides a more robust and clinically reliable solution for automated cardiac MRI orientation labeling, supporting its integration into efficient and scalable clinical workflows.

Introduction The superior longitudinal fasciculus (SLF) is a major association fiber bundle implicated in cognition, visuospatial attention, language, and motor control, and its impairment is linked to several neurological and neuropsychiatric disorders. This proof-of-concept study was performed with three main objectives in healthy adults. First, to fuse whole brain spectroscopic (WBSI) and diffusion MRI (dMRI) derived parametric maps along the SLF I and II segments to quantify their spatial concordance, second, to evaluate regional metabolite concentrations and microstructural properties along these trajectories and finally, to determine the relationships between the WBSI and dMRI parameters within these segments. Methods Ten healthy adults (4F, 6M; mean age 31.4 {+/-} 7.53 years) underwent 3T MRI including multi-shell high angular resolution diffusion imaging and WBSI. After preprocessing and non-linear co-registration, WBSI-derived white matter metabolite maps and neurite orientation dispersion and density imaging (NODDI) / diffusion tensor imaging (DTI) derived parametric maps were spatially aligned and projected along the centroid of reconstructed SLF I and II segments divided into 20 discrete, anatomically contiguous sections. Results A strong spatial alignment between WBSI and dMRI imaging modalities was confirmed by mutual information and Pearson's correlation analyses. Intra-subject repeatability, as assessed from a single participant scanned three times, demonstrated high tract reconstruction reliability (mean Dice similarity coefficients >0.79; track density-weighted Dice >0.97) and acceptable intra-subject coefficients of variation. Inter-subject coefficients of variation were within acceptable ranges ({approx}3-17%) for most parameters, with free water fraction (fiso) exhibiting relatively higher variability. Single and multivariate regression analyses revealed significant associations between WBSI and dMRI tract profiles: choline/creatine (Cho/Cr) and choline/ N-acetyl aspartate (Cho/NAA) ratios showed positive linear associations with intra-cellular volume fraction (ficvf) and fractional anisotropy (FA), and negative associations with mean diffusivity (MD) along bilateral SLF I, with ficvf and MD identified as the strongest combined predictors of metabolite ratios. Conclusion Co-localization/fusion of WBSI and NODDI/DTI data into one framework offers a reliable, user-independent way for mapping regional metabolite and microstructural alterations along the path of SLF. Moving forward, this image processing pipeline has the potential to enhance diagnosis and clinical assessment of neurological disorders linked to SLF damage.

Background: Automated brain volumetry is increasingly used for clinical and research assessment, but its outputs depend on acquisition and reconstruction settings. We evaluated how parallel-imaging acceleration and vendor deep-learning reconstruction (Siemens Deep Resolve) affect automated volumetry on T2-weighted turbo spin-echo MRI, and whether deep-learning reconstruction stabilises measurements across protocols. Methods: Four healthy volunteers each underwent 14 T2-TSE-TRA acquisitions on a 1.5 T scanner, covering seven protocols (baseline, GRAPPA R = 2, 3, 4, and SMS x2,x3, x4), with each protocol acquired separately once with Deep Resolve off and once with Deep Resolve on, for 56 acquisitions in total. Six pipelines were applied: SynthSeg, OpenMAP-T2, and GOUHFI 2.0 for subcortical parcellation, and TotalSegmentator MRI, HD-BET, and NV-Segment-CTMR for whole-brain masking. A three-rater STAPLE consensus of six subcortical structures on the baseline acquisition served as the expert reference. Reference-relative accuracy (absolute percentage error), cross-protocol dispersion (coefficient of variation), and spatial reproducibility (Dice coefficient and 95th-percentile surface distance) were compared between reconstruction states. Results: Deep Resolve improved reference-relative accuracy most for GOUHFI 2.0 (pooled median absolute percentage error 28.4% to 21.9%), modestly for SynthSeg (16.2% to 15.5%), and had a mixed, structure-dependent effect for OpenMAP-T2. It approximately halved cross-protocol dispersion for all three parcellation tools (coefficient of variation 2.62% to 1.39% for SynthSeg, 4.42% to 2.14% for OpenMAP-T2, and 21.75% to 7.87% for GOUHFI) and improved spatial reproducibility in parallel. It eliminated the systematic GOUHFI under-segmentation seen at high acceleration under conventional reconstruction (median deviation -31.6% at GRAPPA R = 4; 14 failure events under conventional reconstruction, none with Deep Resolve). Whole-brain masking was robust in both states. Conclusions: Deep Resolve improved the robustness of automated brain segmentation in accelerated T2-weighted MRI, reducing cross-protocol variability and improving spatial reproducibility, with the largest reduction in reference-relative error in the most acceleration-sensitive pipeline. The elimination of severe segmentation failures observed under conventional reconstruction highlights its potential to support reliable volumetry at higher acceleration. Reconstruction optimization therefore offers a practical route to more consistent quantitative measurements across acquisition protocols.

PurposeWith growing interest in pediatric MRS to understand the neurochemistry of brain development, it is essential to establish pediatric-specific metabolite T1 and T2 relaxation times for reliable quantification of metabolites. Furthermore, characterization of age-related changes in relaxation times may provide insight into developmental differences in tissue microstructure and cellular composition. Therefore, the goals of this study were two-fold: (1) to estimate metabolite T1 and T2 relaxation times in a pediatric cohort; and (2) to assess linear relationships between age and relaxation times. MethodsTI- and TE-series data were acquired from the bilateral thalamus in 20 and 25 participants aged between 0 and 7 years, respectively. Data were processed and modelled in Osprey 2.5.0. Metabolite T1 and T2 relaxation times of total N-acetylaspartate at 2.0 ppm (tNAA2.0), total creatine at 3.0 ppm (tCr3.0), total creatine at 3.9 ppm (tCr3.9), total choline (tCho), myo-inositol (mI) and the combined glutamate and glutamine (Glx) were estimated by fitting the modelled signal amplitudes to TI-series inversion recovery signal and TE-series exponential decay equations, respectively. Linear relationships with age were tested for metabolite T1 and T2 relaxation times using Spearman correlations. ResultstNAA2.0 and Glx showed significant changes in T1 with age. tNAA2.0, tCr3.0, tCho and mI showed significant decrease in T2 with age. ConclusionSeveral metabolites showed significant changes in relaxation times with age in the thalamus region. Therefore, use of age-specific relaxation correction for pediatric studies using the linear equations provided is encouraged for reliable quantification of metabolites.